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0.  Summary 


If  (U,V)  Is  uniformly  distributed  over  the  unit  circle, 

2  2 

then  U  +  V  is  uniformly  distributed  over  [0,1]  and  is 
independent  of  U/V.  This  simple  result  leads  to  an  improved 
method  for  generating  normal  random  variables  via  polar  coordinates, 
and  also  serves  to  relate  normal  probability  measure  and  Lebesque 
measure . 


1.  Main  Paragraph 


Let  (U,V)  be  uniformly  distributed  over  the  unit  circle: 

[  (u ,  v)  !  u^  +  v^  <  l]  •  Then  is  uniformly  distributed  over 

[0,1]  and  is  Independent  of  U/V.  The  proof  is  simple  and  omitted. 
It  follows  that  is  independent  of  [1  +  (v/u)^]"^  and  of 

[1  +  (U/V)^]"^  =  V(U^  +  and  hence 

U  +  V 

(1) 


X  =  v[-  ^  t ]* 

TT^  ■  1T^ 


2  2  -i 

are  independent  standard  normal  random  variables,  since  U(U  +  V  )” 

and  V(U^  +  V^)  are  distributed  as  cos  0  and  sin  0,  and 

[-  2  ln(U^  +  V^)  is  distributed  as  p,  in  the  polar  representation 
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2.  Generating  Normal  Variables 

[1]>  Box  and  Muller  point  out  the  elementary  fact  that  a  pair 
of  normal  random  variables  may  be  produced  in  a  computer  by  generating 
p  and  0,  then  putting  X  =  p  cos  0  and  Y  =  p  sin  0.  The  obvious 
method  -  use  two  independent  tmiform  [0,1]  random  variables 
and  put  X  =  cos(2ttU^)(-  2  In  Y  =  sin(2TTU^)(-  2  In  requires 

cosine,  sine,  logarithm,  and  square  root  subroutines.  This  method  is 
not  very  practical  unless  some  tricks  are  used  to  speed  it  up.  One 
line  of  Improvement  is  to  generate  cos  0  and  sin  0  in  the  form 
U(U^  +  V^)  ^  and  V(U^  +  V^)  as  suggested  by  von  Neumann  [2]. 

The  representation  in  (l)  goes  one  step  further,  taking  advantage  of 
the  fact  that  is  independent  of  U(U^  +  V^)~^  and  of 

2  2  —•A' 

V(U  +  V  )  ^,  is  itself  uniform  [0,l],  and  may  be  used  to  form  p 
by  way  of  p  =  [-  2  ln(U^  +  V^) ]^. 


Although  there  are  faster  methods  for  generating  normal  variables, 
see,  e.g.,  [3]  -  [5],  the  method  suggested  by  relations  (l)  may  be 
suitable  for  situations  where  ease  of  programming  is  the  primary 
consideration.  Furthermore,  a  slight  modification  of  the  procedure 
may  be  used  to  dispose  of  the  problem  of  handling  the  tail  of  the 
normal  distribution  in  one  of  the  super-fast  programs.  Let  r  >  0 
be  a  constant.  Then  putting 


(2) 


Y  =  -  2  ln(U^  +  V^) 


will  produce  a  normal  pair  (X>Y)  conditioned  by  >  r^.  If 

we  want  to  produce  a  nonnal  variable  Z  conditioned  by  |z|  >  r,  we 
may  generate  X,y  according  to  (2)  and  put  Z  =  X  if  |x|  ^  r.  If 
|x|  <  r,  test;  is  |y|  >  r?  If  yes,  put  Z  =  Y,  if  no,  generate  a 
new  pair  X,Y  and  try  again. 

3*  Normal  Measure  "*  Lebesque  Measure 

Finally,  expression  (l)  is  convenient  for  relating  normal 

probability  measure  and  Lebesque  measure  in  the  plane.  Let 

2  2 

be  the  mapping  of  the  unit  circle  C:  [(u,v) :  0  <  u*^  +  v^  <  l}  onto 
the  plane,  R^,  of  points  (x,y)  given  by  the  delations  in  (l) ; 


4 


That  is ,  we  may  find  the  standard  normal  probability  measure  of 
regions  A  by  mapping  A  into  the  unit  circle  according  to  relations 
(2),  then  finding  the  area  of  T(a) . 


This  type  of  transformation  has  been  used  before,  [6].  We 
have  a  simpler  version  here,  which  may  be  put  in  this  form;  To  find 
the  standard  normal  probability  measure  of  a  region  A,  map  each 
polar  coordinate  point  (p,9)  of  A  into  the  polar  coordinate  point 
(e"^  ,0)  of  the  unit  circle,  and  find  —  times  the  area  of  the 
transformed  region.  Inverting  the  mapping  X  =  cos  2ttU^(-  2  In  U^)^, 
Y  =  sin  2ttUj^(-  2  In  1^2^^  gives  a  more  suitable  method  for  converting 
normal  to  Lebesque  measure,  however.  The  procedure  runs  as  follows: 

To  find  the  normal  probability  measure  of  a  region  A,  map  each  polar 
coordinate  point  (p,0)  of  A  into  the  point  (u,v)  of  the  unit 
square  ((u,v) :  0  <  u  <  1,0  <  v  <  l]  according  to  the  relations 


(4) 


u  =  0/2Tr 


then  find  the  area  of  the  transformed  region.  We  have  found  relations 
(4)  to  be  the  most  suitable  for  a  general  purpose  procedure  for  trans¬ 
forming  normal  to  Lebesque  measure,  the  advantage  being  that  it  is 
relatively  easy  to  assign  an  equally  spaced  set  of  points  Uq,Uj^,... 
and  hence  use  one  of  the  standard  numerical  integration  procedures . 
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The  transformation  may  be  described  graphically  in  this  way: 


2^1  277- 


Use  of  Simpson's  rule  gives  the  approximation 

1*(A)  =  MA«)  S  +  4*1  +  2*2  +  4*3  +•••+  2t2„.2  +  4*2„_1 

-^^(2Tru^) 

where  t^  =  e  -  e  ’  ^i  “  “/2tt  +  i<J,  and 

(5  =  (p  -  a)/4vn. 


2n-’ 


This  elementary  procedure  makes  it  possible  to  -write  a  single 
computer  program-  which  will  handle  a  large  number  of  the  commonly 
encountered  regions  in  the  normal  probability  plane  -  polygonal 
regions,  ellipses,  intersections  of  ellipses,  etc. 
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